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Abstract 

Strange nonchaotic attractors (SNA) arise in quasiperiodically driven systems 
in the neighborhood of a saddle node bifurcation whereby a strange attractor 
is replaced by a periodic (torus) attractor. This transition is accompanied by 
Type-I intermittency. The largest nontrivial Lyapunov exponent A is a good 
order-parameter for this route from chaos to SNA to periodic motion: the 
signature is distinctive and unlike that for other routes to SNA. In particular, 
A changes sharply at the SNA to torus transition, as does the distribution of 
finite-time or iV-step Lyapunov exponents, P(Ajv)- 
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Strange nonchaotic attractors (SNAs), which are commonly found in quasiperiodi- 
cally forced systems, are geometrically strange — they are properly described by a fractal 
dimension — but the largest nontrivial Lyapunov exponent A is negative, implying nonchaotic 
dynamics. 

Since they were first described J|, and shown to be generic in quasiperiodically driven 
nonlinear systems || , a number of characteristics of SNAs have been studied, including the 
important question of how they are created. Three main mechanisms or scenarios for the 
creation of SNAs have been advanced. Heagy and Hammel || identified the birth of a SNA 
with the collision between a period-doubled torus and its unstable parent. This mechanism 
is analogous to the attractor-merging crisis that occurs in chaotic systems |4| . Kaneko |5]|| 
has described the "fractalization" of a torus, namely the increasing wrinkling of tori which 
leads to the appearance of a SNA without any interaction with a nearby unstable periodic 
orbit. This route to SNA (and eventually to chaos) has also been observed in several systems. 
Yalginkaya and Lai |7| have recently shown a third route to SNAs via a blowout bifurcation 
||, namely through the loss of transverse stability of a torus. 

We describe a new mechanism for the creation of SNAs. As a function of driving pa- 
rameter, a strange attractor disappears and is eventually replaced by a 1-frequency torus 
through an analogue of the saddle-node bifurcation. In the vicinity of this crisis-like phe- 
nomenon Q the attractor is strange and nonchaotic. We show that the dynamics at this 
transition is intermittent, and the scaling behaviour is characteristic of Type I intermittency 
||. Furthermore, the signature of the transition in A (or in the distribution of finite-time 
Lyapunov exponents) is distinctive and very different from the routes to SNA that have 
been hitherto discussed. The present mechanism is general and is likely to be operative in 
any quasiperiodically driven system. 

For definiteness, we illustrate our results using the quasiperiodically forced logistic 
map || 

x n+ i = a(l + ecos(27T0 n )) x n (1 - x n ) 
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4>n+l = (fin + CO (mod 1). (1) 

where we take u = (v5 — l)/2, the golden mean, and 2 < a < 4. This system is convenient 
to study since the phenomenology is smoothly related to that of the logistic map in the limit 
of e — > 0. Since the driving term is multiplicative, it is clear that the motion will remain 
bound if a(l + ecos(27T0 n )) 6 [0,4], and thus for any a, the largest value of e allowed is 
A/a — 1. It is preferable to work with the rescaled driving parameter, e' = e/ (4/a — 1) since 
this scales the region of interest in parameter space to < e' < 1. 

Fig. 1 is a phase-diagram of the system showing the different possible dynamical 
behaviours — periodic, strange nonchaotic and chaotic attractors, corresponding to the sym- 
bols P, S and C. The different phases are characterised through the largest nonzero Lyapunov 
exponent A = liniTv^oo jf Sili m 1(1 — 2xj)a(l + ecos(27T0j))|, which we calculate typically 



from 10 iterations of the map, as well as the phase-sensitivity exponent ||10|| . The logistic 
map with additive driving || has identical behaviour. 



SNAs occur along the boundary of the chaotic regions |IT], when A takes small negative 
values. In the a-e' plane, the regions of P, S and C behaviour can be interwoven. In 
particular, there are two distinct chaotic phases: the regions marked Ci and C2 which 
are separated by a narrow tongue where the motion is periodic. The behaviour of A as a 
function of a for e' = 1, namely along the upper edge of Fig. 1, is shown in Fig. 2(a). The 
nonmonotonicity of A as a function of a is typical, and several dynamical transitions can 
be described: P — > S, S — > C, and C — > S (note that periodic motion here corresponds to a 
n- frequency torus, with n = 1, 2 or 4; higher n such as 3, 6, 8 etc. can also be observed in 
small windows of parameter space). 

The intermittency route to SNA that we describe in this work occurs along the boundary 
of the chaotic region C2 and the periodic region T on its right, namely at the point marked I 
in Figs. 1 and 2a. An enlarged view is shown in Fig. 2b where the bars indicate the variance 
in the Lyapunov exponent estimated from several computations. There is a transition from a 



chaotic attractor to a SNA at a ^3.405802 with A changing linearly [JES] through zero. The 
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intermittent transition from the SNA to a torus — the phenomenon described here — occurs 
at a = a c = 3.4058088. . .. At this transition, the distinctive signature is an abrupt change 
in the dependence of A on the parameter a, as well as a concurrent marked reduction in the 
variance in A (see Fig. 2b; the bars have been computed from 50 samples of 10 5 steps). The 
disappearance of the SNA or the chaotic attractor is accompanied by intermittent dynamics, 
which can be conveniently studied by co-evolving two trajectories with identical (xq, 0o) and 
e', with different a: since the angular coordinate remains identical, the distance between the 
trajectories is simply the difference in the x n 's. We find that the time between bursts shows 
the scaling r ~ (a c — a)~ e , with 9 ~ 1/2 (the numerical value we obtain is 0.52± 0.04). 
This transition to an intermittent SNA (Fig. 2c) occurs all along the right boundary of C2. 
In the a-e plane, this is a region of nearly constant e. 

The abrupt death of the strange nonchaotic attractor is through a quasiperiodic analogue 
of the saddle-node bifurcation where a period 1 orbit is born. Consider a sequence of maps 
that are rational approximations to Eq. ([[]) by setting uj = uj^ = F^/F^i, where is 
the A;th Fibonacci number. The _Ffc +1 th iterate of the map is a function of x alone (since 
4>F k+1 = 0o)- This allows for the construction of a bifurcation diagram as a function of a for 
fixed e. The case of k = 2, namely u — 1/2 has been studied earlier by Sanju and Varma 
i~3| as the period-2 modulated logistic map. Between the first period doubling bifurcation 



and the 2-band merging crisis, the attractor of the (unmodulated) logistic map consists of 
two branches. (For the logistic map, these bifurcations occur at a = 3 and a ~ 3.678857. . ., 
respectively). The effect of period-2 modulation is to split the two branches of the attractor, 



by shifting one branch to higher values of a ||13|| . The extent to which the two branches 
are shifted relative to one another depends on the value of e. The single attracting orbit of 
period-2 in the unmodulated map becomes two separate attracting orbits for small values 
of the modulation; and depending on how the branches are shifted, namely on the strength 
of the modulation, two very different attractors — for example chaotic and periodic — can 
coexist. For sufficiently large e the two branches can in fact become disjoint, namely the 
periodic orbit of one branch of the attractor is shifted beyond the 2-band merging point (see 



Fig. 2 of Ref. [0). When this happens, the periodic attractor appears as at a saddle-node 
bifurcation, giving rise to a single periodic attracting orbit of period 1. 

It can easily be verified that this mechanism operates with small modifications for the 
higher rational approximations Uk] this is a case of period-Ffc + i modulation JTB). Thus, in 
the limit of — > uj = (\/5 — 1)/2, we suggest that the disappearance of the chaotic attractor 
and the SNA is through this analogue of a saddle-node bifurcation. 

Viewed as a function of decreasing a, this transition shares some of the features of a 
widening crisis in unforced systems [|I| Jl4],|l5| . The discontinuous change of A at a saddle- 
node bifurcation is softened by the quasiperiodic forcing, and for a < a c , the Lyapunov 
exponent shows the scaling (see Fig. 3) 

A - A c ~ (a c - a)", (2) 

with the exponent /i = 0.37 ± 0.03 at e' = 1. To within the quoted error-bar, the same 
exponent is obtained when the probability density in the burst phase is fit to a power law 
|U]-|l6]]. The remnants of the old torus which are now embedded in the SNA constitute 
one of the weakly coupled components of the SNA. These exponents stay approximately 
constant along the C 2 boundary. 

There are at least two other mechanisms for the creation of SNAs that are operative 
in this system. Along the boundary of the chaotic region Ci in Fig. 1, which appears at 
the end of the (truncated) period doubling cascades, SNAs that are formed can appear as 
a result of the collision between a period-doubled torus and its unstable parent [[|. SNAs 
are also created through the process of fractalization || which occurs along the entire left 
edge of the region C2 (for example the cases studied by Nishikawa and Kaneko ||). Along 
the Ci boundary both mechanisms operate and representative points are marked as B and 
A respectively in Fig. 2a. A fractalizing torus gets increasingly wrinkled as parameters are 
varied, and eventually goes chaotic, without any period-doubling bifurcation. These two 
routes can be differentiated from the geometry of the SNA: in the torus collision mechanism 
||, a period 2 fc -torus gives rise to a 2 fc_1 -band SNA, while through fractalization [[|J|, the 
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result is a 2 fc -band SNA. 

The intermittent SNA is similar in many respects to a reentrant phase: the SNA so cre- 
ated is structurally distinct from the SNA born through other mechanisms. This distinction 
also extends to the corresponding tori, and examination of successive iterates of the line 
(1/2, 0) under the map / of Eq. (|TJ) gives an indication of how this occurs. By an extension 



of the e = case [|17]], this line will be attracted to the 1-torus, 



hm/ fc (l/2,0)^T. (3) 

(The attractor of the dynamics, whether periodic, strange nonchaotic or chaotic, is confined 
in a strip in configuration space, the upper boundary of which is the line /(1/2, 0), and the 
lower boundary is the infimum of /■ 7 (l/2,0),j > 2). The convergence behaviour of the line 
to the attractor, quantified through the phase sensitivity index shows that for smaller 



values of nonlinearity and driving, the convergence is slower than for higher nonlinearity. 
Indeed, the torus stabilized at large a, e' (in Fig. 1) is relatively smooth, as opposed to the 
highly structured and fractalizing torus at smaller values of a. 

To further characterise this distinction, we examine the distribution of finite-time Lya- 
punov exponents P(A N ) in Fig. 4, namely the Lyapunov exponent computed from a iV-step 
orbit on the attractor. Although A is negative, for short times the Lyapunov exponent on 



a SNA is positive | lQfl . The distribution P(A^) has not been studied in detail, but it is 
known to have a significant positive tail which does not vanish even for large N. For the 
intermittent SNA, P(A^) decays slowly both as a function of the (local) A^v as well as N. 
This distinction is a consequence of the fact that the variance in the A increases drastically 
through this transition (in Fig. 2b). In contrast, at the other transitions from a torus to a 
SNA (or vice- versa), the variance does not change significantly and with increasing N, the 
P(Ajv) becomes essentially gaussian. The difference in the distributions is not visible for 
N ~ 50 (Fig. 4a) while at N ~ 1000 the distinct character of the intermittent SNA is clearly 
seen (Fig. 4b). These observations are consistent with a coarser interweaving of the positive 
and negative local Lyapunov exponent regions on the intermittent SNA. Furthermore, this 



transition is robust to low-amplitude noise [pL8| , as has been verified by adding a small ran- 
dom component to the x dynamics. The scaling behaviour in A and the probability density 



in the burst phase were observed to hold |L£|] with slightly modified values of the exponent. 

Such quasiperiodically driven systems are interesting from a variety of points of view. A 
potential novel use of SNAs that can be envisaged is for secure communications |20|| : two 
independent SNAs can be easily synchronized since the Lyapunov exponents are negative, 
while the signals from such systems provide effective masking since they appear erratic owing 
to the underlying strange geometry. 

In summary, in the present Letter we have described a new mechanism for the cre- 
ation of strange nonchaotic attractors through intermittency, whereby a chaotic attractor 
is eventually replaced by a quasiperiodic one. In the intermediate region, there are SNAs. 
The distinctive signature of this transition is a sharp change in the Lyapunov exponent 
which shows large fluctuations and scaling behaviour on the SNA side of the transition. 
This mechanism, which is operative at high nonlinearity and large amplitude quasiperiodic 
forcing, proceeds via a stabilization of chaotic motion and is thus like a reentrant phase 
transition. 

The route to intermittent SNA is a general one. We have studied a variety of related 



systems [[HI and find that the phenomenology is very similar: the phase-diagram that is 
obtained for the driven logistic map is typical. Intermittent SNAs always occur in analogous 
regions in parameter space, namely where a tongue of periodic behaviour separates two 
chaotic areas with differing characteristics [^1]]. We expect that in situations where the 
amplitude of quasiperiodic forcing is an conveniently varied parameter for example, this 
mechanism for the creation of SNAs may be experimentally observed. 

This work was supported by grant SPS/MO-5/92 from the Department of Science and 
Technology, India. We thank Sanju for discussions. 
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FIGURES 

FIG. 1. Phase diagram for the forced logistic map (schematic) obtained by calculating the 
Lyapunov exponent as a function of a and the rescaled parameter e' defined as e' = e/(4/a — 1) 
in a 100 x 100 grid. P and C correspond to Periodic (torus) and Chaotic attractors. The shaded 
region along the boundary of P and C corresponds to SNA (marked S). The boundaries separating 
the different regions are convoluted, and regions of SNA and Chaotic attractors are interwoven in 
a complicated manner. The intermittent SNA is found on the edge of the C2 region marked I. 

FIG. 2. (a) The largest Lyapunov exponent A versus a at the upper edge of the phase plane, 
namely at e' = 1 and 2 < a < 4. Ci, C2 correspond to the two chaotic phases identified in 
Fig. 1. T denotes periodic motion (tori); SNAs are found along the boundaries of Ci and C2, in 
the neighbourhood of the points marked A, B and I. (b) An enlarged view of Fig. 2a) near the 
intermittent transition marked I. The bars indicates the variance in A computed from 50 samples 
of length 10 5 . (c) The Intermittent SNA at a = 3.405808 and e'= 1. 

FIG. 3. Scaling of the probability density (ps) in the burst phase (o) and A (□) at the intermit- 
tent transition at I, for a c = 3.4058088 . . . and A c = —0.0283 at e' = 1. The measured exponents 
are « 0.37. 

FIG. 4. (a) Distribution of A^-step or local Lyapunov exponents, computed from N = 50 step 
segments of a long trajectory for the intermittent SNA at the point marked I contrasted with SNAs 
(born through other mechanisms) at the points marked A and B in Fig. 2a. These correspond 
respectively to I (+), A (V) and B (□). (b) Same as (a), except that N=1000. 
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